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We have calculated the real part x' °f the nonlinear dielectric susceptibility of amorphous in- 
sulators in the kHz range, by using the two-level system model and a nonperturbative numerical 
quantum approach. At low temperature T, it is first shown that the standard two-level model 
should lead to a decrease of x' when the measuring field E is raised, since raising E increases the 
population of the upper level and induces Rabi oscillations cancelling the ones induced from the 
ground level. This predicted i5-induced decrease of x 1S a f odds with experiments. However, a 
good agreement with low-frequency experimental nonlinear data is achieved if, in our fully quantum 
simulations, interactions between defects are taken into account by a new relaxation rate whose 
efficiency increases as Ve, as was proposed recently by Burin et al. (Phys. Rev. Lett. 86, 5616 
(2001)). In this approach, the behavior of x' a t 10w T is mainly explained by the efficiency of this 
new relaxation channel. This new relaxation rate could be further tested since it is shown that it 
should lead: i) to a completely new nonlinear behavior for samples whose thickness is ~ 10 nm; 
Oh. ii) to a decrease of nonequilibrium effects when E is increased. 
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Amorphous materials exhibit universal anomalous properties at low temperature. In 1971, Zeller and Pohl JlJ 
\ discovered below 1 K a quasilinear behavior of the specific heat in a number of glasses contrasting with the Debye 
law of crystalline materials. Anderson, Halperin, Varma |2j and Phillips j^] proposed an explaination based upon the 
existence of localized two-level systems (TLS). Their origin may be due to the tunneling of atoms or groups of atoms 
between two equilibrium positions separated by a narrow energy barrier featuring asymmetric two-well potentials. 
They are assumed randomly distributed in energy splittings and tunneling barriers as a consequence of the structural 
disorder of these materials. This model has proven to be successful to understand most salient experimental features, 
^vq ■ The standard TLS model assumes defects do not interact with one another. However, defects are strongly coupled 
to their environment and can emit or absorb phonons. It leads to an indirect interaction between nearest neighbors 
via the phonon field Certain recent failures to explain nonequilibrium data (at a few kHz) |HJ underscore the likely 
involvement of these interactions below 100 mK. However, these nonequilibrium effects are small corrections of the 
kHz stationnary response, and, up to recently, examples of stationnary susceptibilities strongly affected by interactions 
were very rare : in the kHz regime, it was argued that the ultra-low-T (T ~ 1 mK) plateau of the dielectric constant 
in the linear regime, strongly different from the expected logarithmic increase, resulted from interactions ||. Very 
recently, such a conclusion was drawn from internal friction experiments |?]]. 
£j \ In this work, we show that including interactions in the TLS model with a recently proposed mechanism || strongly 
affects the nonlinear stationnary dielectric susceptibility x' of amorphous insulators at a few kHz. A very complete 
set of such data was published a few years ago by Rogge et al. ^ , twenty years after the pioneering work of Frossati 
et al. flipf . In the linear regime, x' decreases when T decreases, reaches its minimum at T rev and then increases 
below T rev (before reaching the above-mentioned ultra-low-T plateau x'piat)- According to the standard TLS model, 
the x' decrease above T rev is due to the progressive freezing of the diagonal (or relaxational) part of the response, 
while the x' increase below T rev comes from the induced off-diagonal (or resonant) part of the susceptibility : this 
effect enlarges as T decreases as do all quantum effects. However, due to the quantum nature of x' below T rev , one 
expects x' to be strongly depressed by a strong measuring electric field E at a given T . This can be guessed from the 
quantum saturation phenomenon which is very general in two level systems |27j. Indeed, increasing E decreases the 
population difference between the two energy levels : as the Rabi oscillations produced by E on the upper level are 
in phase opposition with respect to the ones produced on the ground level, the quantum response, once averaged on 
many independent TLS's, tends to zero when E is increased. Strikingly, Rogge et al. experiments show the opposite 
trend : x'(T < T rev ) increases when E is increased. 

As it is carefully explained in Ref. |J, this behavior does not result from heating of the sample by E. To give 
a supplementary argument with respect to Ref. ||, let us note that if Eu n is the upper field below which the 
dielectric susceptibility is measured as being field independent, one expects that the heating of the sample, for a given 
E > Eii n , is more important when T decreases. A heating effect is thus expected to strech the x'(T) curve of an 
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amount increasing as T decreases, i.e. one expects 



(1) 

E<E Hn 

to hold at low T, i.e. mainly at T < T rev _, As can be seen, e.g. on the Fig. 3 of Ref. ||, the trend of the data is 
exactly the opposite of Eq. (1). Finally, since heating effects can be ruled out, the fact that below T rev , x'{E > Eu n ) 
does not behave as expected from the quantum saturation phenomenon seems extremely intriguing in the framework 
of the standard TLS model. 

However, this was not pointed out since the non linear effects in the TLS model were, up to now, only calculated 
by using the adiabatic approximation jflj. Such an approximation states that TLS's are at every moment at thermal 
equilibrium, i.e., it disregards any coherence effects. It predicts an increase of x' with E, i.e. it qualitatively accounts 
for the experimental behavior. However, in the specific case of the real part of the susceptibility, the consistency of 
the adiabatic approximation is questionnable |l2[ . Indeed, as it is very clearly stated in Ref. JHJ , this approximation 
does not hold for TLS's whose Tunneling energy Ao is too small, and yet it finds that the nonlinear part of x' is 
dominated by the smallest Ao values (see after Eq. (3.30) in Ref. More precisely [|ll|, with po — 1 D the TLS 

dipole, even for the lowest electric fields E ~ 1 kV/m of frequency u) ~ I kHz, the adiabatic approximation fails when 
Ao \/HlupqE ~ 3 /iK, while it is well known, from instationnarity experiments that smaller Tunneling energies 
exist in glasses. Besides, the second puzzling point is that, according to the authors themselves 0|, the reason of 
the increase of x' with E in the adiabatic approximation is physically obscure, which leaves unsolved the question of 
the expected "quantum saturation effect" above mentionned. Finally, several predictions of Ref. JllJ] are somehow 
contradicted by experiments ^ : instead of the predicted T rev oc E 1 with 7 > I, the measured data yield 7 < f/2; 
below T rev , at a given E, the predicted peaked behavior of dx' /dT is not observed; at very low T, the observed E 
dependence of x' p i a t contradicts the predictions. 

This work goes beyond the adiabatic approximation, even though, due to the few simplifying assumptions that 
we have made (see Eq. (2)), we do not intend to yield a fully "from first principle calculation". The key point is 
that phase coherence is not discarded here since non linear effects are treated by a fully quantum non pcrturbative 
method. In the first part, we show that the standard TLS model cannot explain the low-frequency experimental data 
below 100 mK since it yields, at low T, the above-mentioned quantum saturation phenomenon. In a second part, 
interactions between defects are added by using an interaction mechanism proposed very recently by Burin et al. [^) , 
and a successful agreement is obtained with experiments. Finally, we briefly discuss experimental predictions implied 
by Burin et a/.'s interaction mechanism. 
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I. STANDARD TWO-LEVEL SYSTEM MODEL 



A. Bloch equations of TLS 



1. Dynamics of a unique isolated TLS 



Consider a TLS that is sitting in a double- well potential and assume this defect has a dipole moment po . Its energy 
splitting e is related ]l5| to the asymmetry energy A and to the tunneling matrix element Ao, describing transitions 
between the wells, by e = \J A 2 + A 2 ,. Due to finite Ao, the eigen states extend over both sides of the TLS, and the 
position operator r is no longer diagonal in this eigen basis. As a result, when an external electric field E is applied 
to po, the coupling Hamiltonian gE.r is not diagonal in the eigen basis |5| (upon which all the operators of this work 
are expressed), yielding a total Hamiltonian : 



H= \[\l - e ) + ( ^ Ja ) Po-Ecoswt, 

or H = — s.f2, with s = |S where X are the three Pauli matrices and ft is an external effective field (fl components 
are given below, note Q y = 0), which shows an effective spin operator s is associated to the TLS. The systematic 
use of "spin" language comes from the fact that the three Pauli matrices, combined with the identity matrix, form a 
general basis for TLS's. Whatever its physical nature, any operator can be expressed as a linear combination of these 
four matrices, e.g., the density operator p can be written : p = (1/2)1 + (1/7t,)S.S, where S is the quantum mean 
value of the spin operator s . This shows that S x and S y describe the coherence effects contained in the off-diagonal 
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terms of p, while S z is proportional to the population difference between the levels (the occupation probabilities are 
given by the diagonal terms of p). 

The movement of po and thereafter the dielectric response of the material stem from the dynamics of S. For a 
perfectly isolated TLS (note that this implies that T = 0) the evolution of S in the external field S~2 is only a precession 
around the external field J7, as can be seen from the Schrodinger equation which leads fa] to dS/dt = S x S7. 



2. Dynamics of an ensemble of non-isolated TLS's 



At finite T, the dynamics of the TLS must include the relaxation toward its equilibrium value since each TLS 
interacts with its environment (phonons or neighboring defects). Since these interactions occur randomly for a given 
TLS, the dynamical equation must deal with ensemble averaged properties S, i.e. with quantities averaged over many 
similar TLS's. This evolution is given by the Bloch equations, namely 



f =SxO| S <l>reto ; (2) 
C£ Trelax 

where the last term states that the relaxation of S toward the environement equilibrium values < S > re iax must 
be added to the quantum dynamics (see Appendix A). In Eq. (2) it is assumed that the relaxation of a given S 
component, say S x , occurs with a well defined time constant, say t x . In the important case of short time scales, one 
needs to go beyond this approximation since echo signals do not generally decay as a simple exponential (Jl|, 
This subtle effect is irrelevant here since, as already stated e.g. in Ref. jq], we are only interested in the long time 
range solution of Eq. (2), namely x'(lkHz), i.e. we focus on the particular case u>T2 <C 1 (see below). Similarly the 
relaxation term of Eq. (2) might become more complicated in the case of very strong fields |36|, leading, e.g., to a 
S y /T Xj y term in the relaxation of S x (see Appendix B). However this should not be the case here since we only focus 
on the onset of the non linear regime (jiqE will not much exceed k B T). As a result, the relaxation terms can be 
derived quite simply, as we show now. 



i) Phonon induced relaxation. 

Let us first focus on phonon field relaxation. The occupation probabilities are altered by the emission or the 
absorption of phonons, yielding ||l5| a relaxation of S z , with the relaxation time t\ = ki/(cAq) tanh 2 k B T < wnere K i 
is a sample-dependent constant. Since phonon processes occur randomly and independently for various TLS's, they 
break the phase coherence of the ensemble of (noninteracting) TLS's, yielding a relaxation time 2ri for S x and S y . 
What are the thermodynamic values < S Xl y :Z > to which S x>VtZ relax ? By second order expansion of dynamical 
correlation functions, it was shown that this relaxation occurs towards the so-called "instantaneous equilibrium 
values", namely, < S XtVlZ (t) >= Tr(< p(t) > S X} y, z ) where < p{t) >= exp(~H(t)/(k B T))/Tr(exp(~H(t)/(k B T)) is 
the "instantaneous" thermodynamical density operator and k B is Boltzmann's constant. For this result to be valid, 
several conditions must be fullfillcd, among which the most stringent one is, by far : |po.E|r c < Ti where r c is the 
correlation time of the random electrical field acting on a given TLS due to its small interactions with its neighbors 
(see next paragraph ii) ). Finally, these phonon processes yield in the Bloch equations a term (S z (t)— < S z (t) >)/ti 
for the population relaxation, and (S x . y (t)— < S x . y (t) >)/(2ri) for the relaxation of the coherence terms. 

Does Ti depend on time ? On one hand, under the above stated assumption |po-E| t c < h, it was argued [jl6| that 
Ti does not depend on time (see also Ref. fl7| ). On the other hand, one may argue 111] that, since the applied electric 
field modulates the asymetry energy A, one should use 7i(t) = k,\ / (e e / / Aq) tanh ^J S T , where A e ff = A +pqE cosuit 



and e e ff = y ^Iff + ar i se from the diagonalisation of the total t dependent Hamiltonian H. The use of Ti(t) is 

natural within the frame of the adiabatic approximation JlT[ where the system is assumed to be at thermal equilibrium 
at every instant. In our fully quantum approach, the question is much more difficult. In the particular case of the 
low frequency real part of the susceptibility, however, one can easily explain why using either t\ or T\(t) lead to very 
similar results. Indeed, T\(t) and t\ mainly differ only for the TLS's whose gap lie in the range e < pqE. But, as it 
will be shown in the insets of Fig. 1 and Fig. 3, the gaps of the TLS's driven in the nonlinear regime by a given E 
extend on a much larger domain (see sections LB) and II. B)): this is one of the main results of our fully quantum 
approach. Thus the possible time dependence of T\ is not expected to change the results. This was carefully checked by 
performing all the calculations reported here twice, once using n, once using T\(t): the resulting differences between 
both assumptions turned out in any case to be totally negligible. Hence, throughout the paper n is considered as 
time independent, by simplicity. 

With the above relations, we get for the diagonal elements < pi,i(t) > and < P2,2(i) > : 
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< Pi i CO >= — I , tanh v x z -. 

■ y ' 2 2y /Ql + f>2 2k B T 

and 



< 92 2 CO >=- ; tanh z . 

For its off-diagonal elements, it is found : 



where 



r> ^ ^oPo-E 

Siz (t) = —2 — COS Wt, 



o M e o A p .E 

S2r (t = — — — 2 — coswi. 

w h eh 

Finally, one finds for the phonon field contribution: 



,S X >= ™* tanh^±M 
2^ + ^ 2fc B T 



2y/ni + ni ' 2k B r 



and < S y >= 0. 



ii) "Spin-spin" induced relaxation 

Let us now turn to "spin-spin" interactions : for a given TLS, the effects of thermal transitions of its neighboring 
TLS's can be modeled as a small (fluctuating in time) electric field, i.e., as small fluctuating terms SH(t) -C e, fc^T. 
The latter inequality ensures that the relaxation of the population of the levels (involving S z ) will not be sensitive 
to SH(t). It is shown in the Appendix A that, for a given TLS, the oscillations of S X:V (t) are no longer regular but 
progressively deformed by the random 8H{t) terms: due to the absence of correlations between the SH(t) values seen 



by various TLS's, ensemble averaging leads, by cancelation of phases of many TLS's 21 , to a relaxation of S x . y to 
zero (while S x , y remains finite for any given TLS). This happens on a short characteristic time scale T2 <^ O.lw -1 and 
yields a supplementary <S Xj2/ /t2 for the relaxation of the coherence terms. 

The temperature dependence of T2 is not clear at present : in echo experiments ]2S|| , p9| , both T2 oc T -1 as well 
as T2 oc T~ 2 were reported p(J. This might come both from the fact that accounting for the detailed shape of echo 
signals requires a very subtle theory (see e.g. and from the fact that several mechanisms contributes to t%. 

Indeed, the pioneering work p2| of Black et al. predicted a T2 oc T~ 2 dependence but very recent calculations |2^| 
based upon the mechanism used in part II found that T2 oc T _1 could be justified at low T . Since this new mechanism 
will be used in the last section, we use throughout this work r 2 = K2/T, where k 2 is a sample dependent constant. 
In order to try to take into account the various mechanisms which might contributes to T2, the parameter K2 will be 
widely varied, as can be seen in Fig. 2. Last, owing to the smallness of the poE values considered here, we neglect 
any E effect on t-i as explained in Appendix B. 



Hi) Final form of the Bloch equations 
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Inserting the above relaxation terms in Eq.(2), the three Bloch equations can be written as follows: 



n » (t) S y + - _ x + -Z = 0, (3a) 



^-Q x (t)^ + ^(i)5 x + ^ + ^=0, (36) 
at Lt\ T2 

— - + n x (t) s v + — — = o, 3c 

dt n 

where all the S/r terms come from the relaxation processes, while all the VLS terms arise from the quantum 
dynamics, i.e. from the fact that H and s do not commute. 
Equations (3a) and (3b) also write 

^± + zfi z (t) s + + % = in x (t) s z + ^|^, (4) 

at t 2 ZTl 

with 

S+ = S x + iSy 7 

and t* = 2T1 T 2 . 

Let us note that r 2 appears due to the existence in eqs.(3a)-(3b) of the two terms S Xty /{2Ti). Even if they are 
required by consistency (see above and Ref. p4j), these two terms do not exist in the pionneering works accounting 
either for the small instationnarities || or for echo experiments , [^J , ||(| . In fact these two terms play a negligible 
role in the nonlinear susceptibility. To show this, let us first note that as long as T\ > t 2 , one gets r 2 — t 2 , i.e. the 
Eqs. (3a)-(3c) amount to the simpler Bloch equations used before (especially in pulse echo experiments). The key 
point is that, in the (A, Ao) plane, this domain where t\ > T2 is quite large : it is shown in the inset of Fig. 1 and 
in Ref. |^5| that this domain contains, at least, all the TLS's such that e < ei.2 = {k\T '/ ' K2) 1 ^ ■ As shown in the 
inset of Fig. 1, ei,2 — 0.2 K is much larger than the poE values studied in this work. This indicates that the TLS's 
standing out of the t\ > T2 domain should not be affected by E, i.e. they should be in the linear regime (see Ref. 
p6| ). To summarize, nonlinear effects should come mainly from the t\ > ti region where the two terms S Xy y/(2ri) 
are negligible. This will be analytically demonstrated in section B)2J. 

3. Non perturbative resolution of the Bloch equations 

The Bloch equations cannot be solved analytically and even their numerical resolution is so far a great challenge. 
However, in the audio-frequency range, some approximations can be made which strongly simplify the calculations. 
As r 2 is much shorter than the typical time (~ ^r) to modify the populations, S z may be considered constant |27]] 
in the right hand-side of Eq. (4). The coherence terms follow adiabatically the population evolution. They reach at 
every moment the stationary state corresponding to the "frozen" occupation numbers. 

Therefore, Eq. (4) can be solved independently of Eq. (3). The stationary solution of Eq. (4) is 

- _ in x S z + <S X > /2ti 

in, + 1/75 ' () 

which inserted into Eq. (3) leads to a differential equation for S z : 



dS z n§/rg - s z ~<s z > Q X Q Z <s x > 
dt + ni + i/r* 2 z+ n ni + i/r* 2 2 n ' 



S z (t) in Eq. (6) is expanded into its Fourier series to get its stationary state. The expansion is limited to a finite 
number of harmonics. This number, of the order of 10, is found a posteriori when a stable and accurate result is 
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obtained. So the differential equation is equivalent to a linear system whose solutions are the harmonics S™. The 
inverse Fourier-transform gives the periodic evolution of S z (t) . The coherence terms S x and S y are deduced from Eq. 
(5) where S z (t), the solution of Eq. (6), is inserted. Finally, the first harmonics of S x (t) is sought, to be included 
into the dielectric susceptibility (see Eq. (7) below). 
Indeed, the susceptibility [|] of a single TLS reads 

X=4^cos^f + ^|Y (7) 



E| \e h e h , 

and it must be averaged over the distribution of TLS's || and over the dipole-orientation angle 9 to yield the total 
susceptibility of the sample : 



X = P dA — ° / d(co S 0)x(A,A o ,0). (8) 

JO JA 0ml „ A J-l 

In the remainder of this article, we concentrate on the real part % of x which is linked to the capacitance of the 
sample, i.e., to its dielectric constant e r by : 

i *' 
e r — 1 = — . 



B. The quantum saturation effect: the quantum part of x'(T) is depressed by a E increase 



1. Numerical results 



We have used the standard values for amorphous-SiC^: po = 1 D, P = 3 x 10 44 Jm -3 , k\ = 10 -8 sK 3 (all the 
energies in n taken in K), A 0min — 10~ 6 K, A 0max = 10 K, A max — 10 K. As explained above, we took r 2 = K2/T, 
where K2 was ranged from 3.10 -11 sK to 10 -7 sK, allowing to check our fundamental assumption UT2 <C 1 provided 
T > 0.5 mK. Last, the numerical relative accuracy of our simulations was, in any case, better than 10~ 3 : this was 
checked very carefully, both by increasing the number of harmonics when solving Eq. (6) and by letting the successive 
integration procedures converge to better than 10 -4 . For each set of parameters E,T, Ki,K2, Ao ;m i„ at least 4 x 10 4 
couples of (A, A ) were computed. 

The simulations are displayed on Fig. 1. The resonant response (low temperature) is strongly depressed by the 
drive level, while the relaxation contribution (high temperature) is little affected. This is at odds with the experiments 
H where increasing E leads to an increase of both the resonant response and of its slope \de' r /dT\ below T rev . Let 
us note that the curve labeled "linear response" was obtained independently by a standard series expansion of the 
Bloch equations keeping only, as in Ref. |5|, the terms proportionnal to E : as E is made very small, the nonlinear 
calculations very precisely converge towards the linear regime. 

However, the extreme sensitiveness of the resonance to the external field is very striking. It decreases rapidly 
while |po-E| << ksT. The low-temperature phase-coherent upturn is destroyed by its environment (the external 
field) , although the perturbation is much smaller than any thermodynamical quantity, which suggests that this effect 
has a quantum origin. This is further confirmed by the inset of Fig. 2 showing the influence of T and t 2 on 
Sx'(E,T) = 1 — x' (E,T)/x'{0,T) : for a given E, the smaller T, the larger 6x', which is expected since quantum 
effects generally increase as T decreases. Similarly, 5x' is larger when K2 is made smaller, i.e., when quantum coherence 
is made more "fragile". Finally, the dimensionless Sx' appears to depend not only on E,T,k,2 but also on Ki, and it 
is shown in the main part of Fig. 2 that all these dependencies are a universal function of a dimensionless scale rj, 
namely, : 



, J 0.1 x Jrj if n < 1 , PoE m 
S * = { 0.1 X ln( ?7 ) if r, » 1 Wlth " = W ' (9) 

where a ~ 0.45 ± .05 and ln(?7) might be replaced by a power law of rj with an exponent lower than 0.1. This 
universal Sx'i'u) dependence holds only when the relaxational part of x' can be totally neglected, i.e., well below 
T rev ~ 50 mK : in Fig. 2, only data corresponding to T < 10 mK have been plotted. For these low T, Sx'iv) remains 
universal even when (ki, K2, E) are varied over several decades. The factor K\ / (T 2 K2) in 77 becomes very large at low 
T, yielding nonlinear effects even for very small E: this expresses that the lower T, the smaller the onset field of the 
nonlinear regime, as already seen on Fig. 1. Let us mention that the data of Fig. 2 correspond to the particular case 
6> = 0. 
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FIG. 1: Inset: At T = 10 mK, Ki = 10~ 8 sK 3 , and k 2 = 10~ 8 sK, the domain of TLS's such that r 2 < Ti is quite large and 
contains all the gaps smaller than ei,2 = (kiT/k2) 1/ ^ 3 -see Even at po-E = 5.12mK this domain is larger than the one of 

the TLS's driven in the nonlinear regime defined by e < e on set — 70 mK (see Eq. (11c) ). Note that e on set 3> PoE (poE is the 
small black area very near the origin) : this explains that the nonlinear effects are visible even at very low fields, as shown in the 
main Figure. Main Figure: Dielectric susceptibility of amorphous- SiO 2 at 1 kHz vs temperature simulated at various fields 
-the value of poE in Kelvin labels each curve- within the standard two-level system model with the following set of parameters: 
po = lD,m = 10" 8 sK 3 , k 2 = 10" 9 sK, A , mi „ = 10" 6 K, A max = A 0mal = 10 K, P = 3 x 10 44 Jm -3 . The low-temperature 
response vanishes rapidly as the electric field is increased due to the quantum saturation phenomenon. The linear response was 
obtained by an independent perturbative method. 
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FIG. 2: Inset : 5x' = 1 - x'( E , T)/x'(0, T) plotted versus pqE (in kelvins). Curve A corresponds to p = 1 D, Ki = 10" 8 sK 3 , 
k 2 = 10" 9 sK, A ,min = 10" 6 K, A ma x = A 0ma x = 10 K, P = 3 x 10 44 Jm -3 and T — 2 mK. The other three curves show 
the effect upon quantum saturation of the parameter which was changed with respect to A : increasing T, as well as decreasing 
Ki, decreases 5x't while decreasing K2 increases 5x\ as expected due to the quantum nature of 8x ■ Main figure : The various 
influences of the simulation parameters can be reduced to a universal function of the dimensionless variable n = JX 1 ) a 

with a = 0.45 ± .05 numerically. The dashed line shows that 6x' oc yjrj when r) < 1. The various parameters were ranged over 
several decades : 10" 10 sK 3 < m < 10" 8 sK 3 ; 3 x 10" 11 sK< k 2 < 10" 7 sK; 10~ 6 K< A , mi „ < 10" 4 K; 10" 8 K< p E < 3 
mK. The data of this figure correspond to the particular case 6 — 0. 
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2. Physical interpretation 

To further understand the universal 5x'(v) an d demonstrate its quantum origin, let us briefly go into the structure 
of the Bloch equations. By using the identity Q x < S z >— fl z < S x >, Eq. (6) can be written : 

dS z S z _ < S z :- - J 77j = n ( 1 + hi+(nl 2 r]y) . . 

* T -~ T - 1 * - *(i + ihSs&) ' ( 1 

In Eq. (10), one gets at E — > : t z = t Zj i = n. As argued in section T)A)2), the nonlinear behavior should come 
from the TLS's such that n > r 2 : in this case we see indeed from Eq. (10) that increasing E decreases r z much more 
than t Zi \. This strongly depresses the off diagonal suceptibility, as we show now. 

Let us first derive, from Eq. (10), the critical value E* such that 1/t z becomes larger than l/r^i: focusing on 
the gaps e lying within the T\ > t 2 domain, i.e. in the domain where r 2 — t 2 , E* is determined by the condition 
TiT 2 r^ ~ l + 02 T | ; y ie idi n g : 

= g R , T if fcser 2 < ?i (11a) 
= ey^ if fc B er 2 > h (lib) ' 

where all the energies are expressed in kelvins. With the standard values K\ = 10~ 8 sK 3 and k 2 = 10~ 9 sK, we see 
that pqE* is much smaller than e. Indeed, for T = 10 mK we get poE* /e = 2.10 -5 for the smallest gaps following 
Eq. (11a), and, for example, poE*/e < 3.10~ 3 for the gaps e ~ k B T which follow Eq. (lib). Solving Eq. (lib) with 
respect to e, for a given E, leads to a characteristic gap 

/ \ 1/4 

Conset = VPoE I , (11c) 

where all the energies are in Kelvins. For the highest pqE ~ 5.12 mK studied here, we get e onS et — 70 mK. As 
shown in the inset of Fig. 1, e onse t is both much larger than pqE and corresponds to a domain smaller than the one 
defined by our assumption n > t 2 . 

To show that E* in Eq. (lib) is indeed the critical field for a given TLS, at which the kind of nonlinearities of 
Figs. 1-2 onsets, let us now compare x'(E -c E*) and x'(E*). 

i) If E <C E*, we get from Eq. (10) t z ~ r Zj i ~ n. Solving Eq. (10) is straightforward and leads for the n th 
harmonics of S z (t) : 

< S n > 

^ = 1 + nW^ (12) 

where < S z > is the n th harmonics of < S z (t) >. Remembering that the region of interest is e < e onset , it 
can be checked that lot\ 3> 1 for basically all the considered TLS's. This yields, from Eq. (12), S z (t) ~< S z >. 
Furthermore, since poE <C e due to Eqs. (11), we get < S z (t) >~< S z >, which, once combined with the identity 
&x < S z (t) >= fl z < S x (t) >, yields S z (t) ~ fi z < S x (t) > /9, x . Once reported into Eq. (5), this yields : 

S.W ^f'ff f fffrg), ( 13 ) 

1 i- S 2 Z T 2 

where in the last factor the fact that fi z r| ^> r 2 /(2ri), which holds for any reasonable set of (k\, k 2 ), was used to 
drop the term r 2 /(2ri). 

ii) For E = E* , we get from Eq. (10), r z .i ~ t\ and Ti/2 < T z (t) < t\. The fact that t z is now smaller than r Zi i is 
responsible for the onset of nonlinear effects. This can be seen by setting t z = n/2 throughout the electrical period. 
With this simplification, one gets, with a derivation similar to the one yielding Eq. (13) : 
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The off-diagonal part of the response in phase with E is x' x ^ &1./E : it is read directly from Eqs. (13)-(14), 
remembering that < S x >oc Ecosut. This yields x' x (E — E*) ~ ^x' x (E -C _E*), where the factor 1/2 comes from 
the above relation r z = ^t z ,i, which was a simplification of the case E = E*. The comparison of Eqs. (13)-(14) 
is thus only semi-quantitative, but it yields the main two features of the quantum saturation phenomenon : first 
x'(E*) < x!{E <C E*), second this effect comes from the off-diagonal part of the susceptibility, i.e., it is purely 
quantum (the diagonal susceptibility x z cx «5 Z /.E is much smaller than x' x due to the fact that ujt\ 3> 1 below T rev ). 

We have here an example of quantum decoherence |53|. It is not surprising that these effects were missed by the 
adiabatic approximation mentionned in the introduction since, in this approach, T2 has disappeared, yielding for the 
nonlinear onset || no other possibility than |po-E| ~ feT, as expected for a system at equilibrium. Moreover we 
have shown that the quantum saturation depends on the precise coupling of the three Bloch equations, i.e. of the fact 
that t z evolves faster with E than r z \ : this is out of reach for the adiabatic approximation since it contains only 
one differential equation Jll| instead of Eqs. (3a)- (3c). Finally, the results of Figs 1-2 do not depend on the precise 
microscopic mechanism involved in T2, but only on the fact, well established by echo experiments, that, for a vast 
subclass of TLS's one has t 2 <ti : this is the main reason of the ^-induced depression of x' of Figs 1-2. 

3. Effect of the density of states 

More can be learnt from Eqs. (11), and more precisely from Eq. (lib) which holds for the vast majority of the 
TLS's responsible for the nonlinear behavior. First, let us note that the onset field E* increases as \/k?Jk\ : this 
suggests that the depression of x\ when E is increased, depends on E^J k±/k2, which, remembering that n\/ K2 is the 
square of a temperature, leads to the dimensionless scale pqE / '(fc^T) <J k± / '(T 2 k%) as the natural parameter for the 
quantum saturation phenomenon. This dimensionless scale matches exactly the definition of 77 in Eq. (9). 

Second, from the above discussion of Eqs. (13)-(14), the TLS's such that e < e onset are already in the saturation 
regime, while the gaps larger than e onse t are hardly altered by E. It is thus natural to consider the number of TLS's 
such that e < e onse t as an estimate of the amplitude of the quantum saturation phenomenon 1 — x'(E,T)/x'(0,T), 
stating : 

1 - X '(E, T)/ X '(0, T) cx £°" S " P(e)de cx sfE cx y/rj, (15) 

where the last equality was obtained by using the above-stated relationship E oc 77; while the second equality uses 
both Eq. (11c) and the fact that the energetic density of states P(e) is a constant due to the standard distribution 
P(A, Ao) = P/Aq. Equation (15) yields exactly Eq. (9) derived from the numerical simulations. This argument 
enables to state that the small corrections to the standard P/Aq previously proposed only yield small changes to the 
behavior of Figs. 1-2 : this is true, e.g., for the slight depression of the density of states at small gaps derived by 
Burin pl[ , as well as for P j 'A +v with \y\ <C 1 proposed in Ref. p2] j. 

To summarize this section I), solving the Bloch equations leads to the quantum saturation effect, i.e., to a strong 
decrease of the off-diagonal part of x' when E is raised. This effect holds for a very large set of k% and n 2 -the 
main parameters of the model-, and it mainly comes from the TLS's such that e < e onset < For an ensemble of 
TLS's with a P/Ao density of states, quantum saturation goes as E 5 , and such an exponent justifies a posteriori 
the nonperturbative character of the method used here. Last, the quantum saturation phenomenon onsets for fields 
E* <C ksT/po, as seen from Eq. (9). It is thus non-negligible since the field is, in most experiments, decreased well 
below ksT/pQ. However, in the literature, the trend of the data is systematically the opposite of the one of Figs. 1-2. 
Since -see Appendix B- more general Bloch equations, corresponding to larger E, should not qualitatively change the 
results of Figs. 1-2, we conclude that the standard TLS model cannot account for the basic features of the nonlinear 
experimental data in the kHz range. 

II. ADDING INTERACTIONS 

A. Burin et al's mechanism 

At this step, at least one drive-dependent parameter must be added into the model to explain the large discrepancy 
with the experimental data. Moreover, it must enhance the relaxation process at low temperature, since coherence is 
broken by the external field as shown in Figs. 1-2. 
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Recently, Burin et al. |B[ proposed an additional field-induced relaxation mechanism. They show that the resonant 
dipolc-dipole coupling, which is so small in glasses, can be strongly increased by a low-frequency electric field. Indeed, 
thermal excitations, which are at zero-field localized on each TLS, tend to delocalize by hopping to resonant nearest 
neighbors. This is due to the fact that resonant hopping demands both TLS's to have very close values of both A and 
Ao : as the electrical field modulates the TLS parameter A, the probability of finding, for a given TLS, a resonant 
TLS, increases from a negligible value at very low E, to a non-negligible value above a threshold of the external field. 
Let us note that this mechanism transports energy; hence it can be treated as a new relaxation mode. 

The frequency must be small for the electric field to have time to modulate the coupling parameters. This is of no 
consequence here, since our crucial assumption uit-i <C 1, leading to Eq. (5), already restricts our work to the low 
frequency case. Another assumption is that the external field amplitude is smaller than the characteristic splitting 
energy ~ kgT, in order to treat the field as a weak perturbation. The typical values of the frequency and |po-E| are 
respectively 100 Hz and 1 mK but may be softened as a rigourous determination is out of reach. 

When the electric field increases, so does the probability of finding a resonant neighbor close enough to yield 
tunneling with not too small a probability : the one-particle excitation will relax more rapidly at high E towards 
another site. One can show the relaxation rate is proportionnal to the square root of the drive level To include 
this new enregy relaxation channel, we set in Eqs.(3a)-(3c) r-f 1 = r^^ h + where T\. v h is the phonon field induced 
relaxation mechanism used throughout section I) and where 



TIB = -7^=, (16) 

V|Po-E| 

with the constant B = 10~ 5 sK 1 / 2 for physically reasonnable parameters ||. As a result, increasing E at any given 
T leads to an increase of the susceptibility \' '■ this shows that Burin et a/.'s mechanism is strong enough to overcome 
the decrease due to the "quantum saturation phenomenon". However, the agreement between the set of calculated 
curves (unreported) and the data is very poor since the net increase of x'(T) when E is increased is stronger at high 
T than at low T. This is due to the fact that, since relaxation dominates the total response, the most influent TLS's 
are such that e < fcgT : their number enlarges with T and so does their supplementary relaxational response due to 
the new relaxation channel t±b- 

To interpolate between Fig. 1 and Eq. (16) which appear as extreme cases, one mig ht use the very general argument 
that interaction effects should disappear at high T, e.g., above 100 mK -see Ref. P0|| -. This demands that the chosen 
t ib(T) becomes infinite, i.e., negligible, at high T. Such a general requirement can be of course modeled by different 
laws but all the ones we tried gave the same kind of behavior for the susceptibility. This is why we report on the 
calculations using a simple law, namely, 



tib (T) = - -t b /t witl1 Tb = 15 mK ' ( 17 ) 

where t\b is given by Eq. (16) and the thermally activated behavior models a dipole-dipole coupling constant of 
Tb = 15 mK : the energy scale Tg can be deduced from Fig. 3 of Rogge et al.'s data || on a-SiO x since x' becomes 
T-independent below 15 mK even for E values ten times larger than the range of the linear regime. Of course, this 
Tb scale can be adjusted empirically since the T where \' becomes T-independent depends on the material. As 
the coupling constant goes as g/ |r — r'| 3 and as ||, for a-Si02, g ~ 10 Knm 3 , we get a mean distance Xb between 
interacting dipoles of nearly 10 nm. 



B. Numerical results 



The modified-model predictions using Eq. (17) are displayed on Fig. 3. The values of |po-E| have been limited to 
10 mK because of the restrictions on both the Bloch equations and the field-induced mechanism. A trend completely 
different from the one of Fig. 1 is obtained at low temperature since an increase of the response is observed when the 
drive level increases. 

By computing separately (unreported) in Eq. (8) the two terms of the right hand side of Eq. (7), we checked that 
x' x behaves qualitatively as in section I) and that the new trend of Fig. 3 is due to the diagonal part x' z - To explain 
this new behavior, one first note that t\b(T) is now the upper bound of n, even for the numerous TLS's whose small 
A value lead, in section I), to a very large n. With ujt\b{T < Tb) % 1, the 1/(oj 2 t 2 ) cutoff of S z seen on Eq. (12) 
has now disappeared, i.e. the dS z /dt term in Eq. (6) can be dropped, yielding : 
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FIG. 3: Main Figure : Simulation of a-Si02 susceptibility at 1 kHz vs temperature with Eq. (8) and the same parameters as in 
the main part of Fig. 1. The calculations were done within a modified TLS model where excitations are no longer localized but 
can experience field-induced hops to neighboring sites, which is modeled by an additive relaxation channel (see the definition 
of tis(T) in Eq. (17)). The data show a linear behavior at low enough drive levels (the poE values label the curves), an 
evolution of T rev with E compatible with experiments and a substantial decrease of the T dependence of x' a t l° w T ; Inset 
: For poE <C ksT, in the (A, Ao) plane, < S\ > is not negligible only within the e < 2T domain. Even for poE = 0.8 mK, 
the hatched area where r z < t z ,i has a non negligible size with respect to this e < 2T domain : this yields a supplementary 
T-dependent contribution to the diagonal susceptibility x'z which overcomes the ^-induced depression of x'x seen on Fig. 1, 
and yields the i5-enhanced \' trend seen on the main part of the figure. 



S z {t)^—{~S z (t)), (18) 

T*,l 

where t z ,t Zi i are defined in Eq. (10). At E — > 0, one has r 2 ~ r z .\ ~ n, yielding with Eq. (18), S z (t) ~ (S z (t)). 
With the additionnal remark that (S*(e < 2T)) ~ hp E / '(Ak B T) while (5^(e > 2T)) ~ 0, one gets, with the standard 
P/Ao density of states, that ~x! z (T) oc +lnT : this is the trend seen above T rev . 

To explain the behavior below T rev , the key point is that for quite a large domain in the (A, Ao) one has t z /t Z) i < 1 : 
since this factor is T dependent, it will modify the T dependence just above derived for x' z horn Eq. (18). Focusing on 
the e < 2T gaps, we get, below T rev , n ~ ti„b(T) and t\ 3> T2 : the condition t z /t z< i < 1 amounts to TiT 2 f^ > O^rf, 
i.e. : 



£ < 2poE t\b{T)/t2 sin^ + cos^J with <j> = arctan-^-, (19) 

The t z /t z ,i < 1 condition is shown, as a hatched domain, in the inset of Fig. 3. Even for the lowest E studied 
here, it is not negligible with respect to the e < 2T area. Since in the hatched domain one has t z /t z i ~ T2& z / '(Tif2|) , 
this factor remains T dependent even below Tb when t\b{T) has reached its maximum value : this is due to the 
fact that T2 remains T dependent even at very low T. With (Sl(e < 2T)) ~ hp E/(4k B T), integration of Eq. (18) 
within the hatched area yields a contribution Sx' z oc E 3 / 4 /T 1 / 2 . Thus: i) this term increases as T decreases; ii) Sy' z 
increases with E, i.e. it can overcome the iJ-induced depression of x'x- Disregarding the slight difference -see |35| - 
between the Sx' x °c —E 1 ^ 2 seen for the quantum saturation phenomenon and the Sx' z oc +i? 3 / 4 , the linear regime 
of Fig. 3, up to poE — 0.32mK can be seen as resulting from the compensation of both effects. At higher E, the 
Sx' z increase dominates over the .E-induced depression of x' x > yielding a net increase of x' w hh E. Note that x'z(E) 
becomes T independent when T < poE/ks ■ in this case, indeed, (Sl(e < 2T)) is no longer T dependent. This yields 
the substantial decrease of the T dependence of x' seen for the two highest E values on Fig. 3. 

Last, the off-diagonal susceptibility x' x °^ $x mainly behaves as in section I), i.e. we recover the quantum saturation 
phenomenon yielding, when E is raised, both a decrease of x'x an d of the slope \dx' x /dT\. With respect to section 
I) the quantum saturation effect is somehow weakened, which can be understood since, for a given E, the number of 
TLS's lying within the e < e on set domain of Fig. 1 is larger than the corresponding one in the hatched area of Fig. 
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3. Finally, the variations of x' x with T remain smaller than the ones of x' z > excepted in the case where T < pqE/Ub '■ 
the small T dependence of x'(T < Tg;p E > 5 mK) is thus the only case where x' x dominates the T behavior of x' 
on Fig. 3. 

To summarize, the bigger the electric field, the smaller the field-induced relaxation time (see Eqs. (16)-(17)), which 
enhances the relaxationnal part of the response, leading to a net increase of x' with E at a given T. At a given E, 
when T decreases below Tb, the x'fT) increase is due to the fact that T2 is still T dependent : this is, of course, out 
of reach for the adiabatic approximation where Ti has disappeared. Finally, inserting Burin et aVs new relaxation 
rate in Bloch equations allow to account for the main trend of the nonlinear data : however, in this approach, the 
so-called "resonant" regime below T rev is not a coherent one but, mainly, a field-enhanced relaxation regime. 

C. Comparison with experiments 

On Fig. 3, one observes a pseudolinear regime up to pqE ~ 0.05/cbT where the dielectric response is quasi- 
independent on the external field. This value of the electrical field agrees with the experimental linear regime, which, 
depending on the materials, extends up to PqE / '(fcgT) in the range [0.02; 0.12] (see Figs. 3-5 of Ref. ||). We checked 
that this pseudolinear regime comes from the form of n t B <x E~@ where (3 takes the highly nontrivial value 1/2. 
Setting lower values for /3, such as (3 = 0.1, yields the quantum saturation phenomenon to dominate, leading to the 
same trends as in Fig. 1, at odds with experiments. Setting (3=1 leads to the tendency of Fig. 3 but with a linear 
regime reduced to poE/(ksT) < 0.01. The second key point is the trend of the reversion temperature T rev with E : 
using Eq. (17), i.e., (3 — 1/2, leads T rev to increase by a factor three when E = 30 x E rev , where E rev is the electrical 
field such that the nonlinearities onset at T rev . This is in good agreement with Fig. 3 of Ref. |l . On the contrary, 
using (3 = 1 leads T rev to increase much faster with E : T rev (E = 30 x E rev ) = 30 x T rev (E = 0). Finally, the key 
role of /3 = 1/2 is somehow reminiscent of Eq. (9) where Sx' °c y/E, even if an analytical argument supporting this 
idea is still lacking. 

With respect to experimental data, a failure, at this step of the discussion, is the ratio between the two slopes 
de' r /d\nT below and above the reversion temperature. In Fig. 3 this ratio is near -1.7:1 instead of -1:1 in most 
experiments. Furthermore, the low-temperature /wear-susceptibility data tend to a T-independent plateau while 
they do not in our simulations. At very low temperature, interactions are likely to be so strong that the independent 
TLS model does not apply anymore, even with a renormalized relaxation time such as that of Eq. (17). A transition 
toward a dipole-glass was invoked to explain the behavior of the samples whose x' no longer depends on T below a 
few mK. In this picture, dipole orientation is progressively frozen, which would lead to a plateau of the susceptibility 
p|, p9|| : by continuity, this would weaken the slope ratio near -1:1. Since the TLS model should not apply at very low 
T, it is not surprising that the plateau of the susceptibility measured in the nonlinear regime is not well accounted 
for by Fig. 3. Indeed, Fig. 3 does not show a completely T-independent plateau but only a substantial reduction of 
the T-dependence of x' at low T: as stated in II)B), this is due to x'x which still exhibits a small T dependence, even 
when x' z nas turned into its T independent regime. However, if, on Fig. 3, the susceptibility is frozen below a given 
T, one gets plateaus for x' whose heights depend on E, as in experiments. Finally, pushing [3 toward 1 strengthens 
the tendency of x' to become T independent at low T (unreported), even if j3 ~ 1 leads to the above-mentioned 
discrepencies with respect to experimental data. Let us note that some materials (see Rogge et al. ^]) do not yield 
any sign of such a glass transition even at T = 0.6 mK. 

D. New predictions 

Let us move briefly to the physical predictions implied by Burin et a/.'s mechanism. Remembering that the inequality 
u>T2 <C 1 allowed the key simplification for the derivation of x'{E,T) -see Eq. (5)- we restrict ourselves to the kHz 
range where this condition is fulfilled. Two main predictions can be done : 

i) t ib(T) will be suppressed in samples whose thickness h is smaller than the distance Xb separating the quasi- 
similar TLS's required by Burin et a/.'s mechanism. Indeed, at distances larger than h, dipolar interactions within 
the dielectric will be suppressed by the screening effect of the numerous electrons of the electrodes. Thus, if h ^ Xb, 
one should observe a non linear behavior such as the one calculated in section I) -see Fig. 1-, where the quantum 
saturation of the levels only remains. In other words, ranging h from a fraction of As to a few As in a series of samples 
and studying x'(E, T) should lead to a gradual transition from Fig. 1 to Fig. 3 if Burin et aVs mechanism is relevant, 
while it should not affect the non linear behavior in the standard TLS model. Note that such an experiment looks 
feasible due to the quite large value of Xb — 10 nm, -sec II)A)-. This is due to the fact that Burin et al.'s mechanism 
requires the two interacting TLS's to have both very close values of A and very close values of Ao : these conditions 
are stringent enough to make As much larger than the distance between a given TLS and its nearest neighbor. 
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ii) The net relaxation frequency r-f + of a given TLS increases as E increases. Thus, nonequilibrium data 
should be of smaller amplitude when E is raised. Indeed, they are currently interpreted as resulting from the very 
large t% existing in any glass due to the subclass of TLS's whose energy barrier is so high that Ao is very small. These 
very "slow" TLS's have an extremely delayed response to any change of the external constraints, such as the d.c. 
electrical, or strain, field imposed to the sample: these TLS's yield an excess of states at low energy with respect to 
the equilibrium density of states, the latter having a small depression at low energies due to TLS-TLS interactions. 
To our knowledge, the influence of E on nonequilibrium phenomena has been reported only once, in Rogge et aZ.'s 
work devoted to nonequilibrium phenomena on a mylar sample |ffl| . Applying a relative strain field T to the sample 
leads to a sudden jump of the dielectric capacity C, measured at 5 kHz, followed by a logarithmic relaxation. At 
T = 11 mK, i.e., well below T rev , and with T = 2.7 x 1CT 6 , the initial relative jump is dC/C = 13 x 10~ 7 if the 
measuring field is E = 5 x 10 4 V/m (see Fig. 1 of Ref. pl|), while it decreases to dC/C — 4.5 x 10~ 7 if the measuring 
field is E = 8.5 x 10 4 V/m (see Fig. 2 of Ref. [pl|). Let us note that, with po = 1 D and a relative dielectric constant 
of 5, E — 5 x 10 4 V/m amounts to an energy of 10 mK, of the order of T : in terms of our Fig. 3 this means that 
one stands just above the pseudolinear regime, i.e., in a regime where our calculations, as well as Burin's mechanism, 
should apply. Even if this was not investigated systematically, this single experimental datum favors the idea that 
nonequilibrium effects should be of smaller amplitude when E is increased, due to the interaction-induced reduction 
of the diagonal relaxation time. 



III. CONCLUSIONS 



In conclusion, we have simulated the nonlinear dielectric susceptibility of amorphous materials by using the TLS 
model. Phase coherence effects have been taken into account, which is the main difference with the adiabatic approx- 
imation. In the kHz range, the standard TLS model yields a nonlinear behavior at odds with experiments due to the 
field induced depression of the quantum response. However, it was possible to fit in many details the experimental 
low-temperature field-induced rising response by adding a new relaxation mechanism based upon the existence of 
interactions below 100 mK. In this approach, the low temperature response mainly loses its quantum origin at low 
frequency. Our work stresses the necessity to inject interactions into the TLS model to get satisfactory predictions. 
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IV. APPENDIX 



A. Phase decoherence induced by small TLS interactions 



In this Appendix, we aim at giving some physical insight into the relaxation term introduced in the dynamics of 
an ensemble of TLS's due to their small mutual interactions. Expanding on the assumption that these interactions 
are much smaller than the other relevant energy scales (such as T or the gap e), the basic idea |ll| is to model these 
interactions by a small random electric field acting on each TLS. This idea is not new |22[ , Jl6[ , and numerical results 
are presented here only to help understand theoretical results. 



1. Interactions effects when the measuring field E = 



Consider first the case where the measuring field E = 0. Modeling mutual interactions between TLS's by a random 
electric field leads, for a given TLS, to a total Hamiltonian given, by : 



H =2 



e 
-e 



A 

in 
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where the electric field ~E ran d is random in time for the considered TLS, and, at a given instant t, varies randomly 
for various TLS's. Note that Eq. (Al) is expressed in the eigen basis of the TLS. 
Defining the density operator p(t) by : 

P (t)=( L 2 +z * + (A2) 
\x - ty ^ - z J 

it is clear that x,y,z are, respectively, the quantum mean values of the three spin operators (S x ,S y ,S z are the 
corresponding symbols once the ensemble average over many similar TLS's is made). By using ihp = Hp— pH, where 
the dot stands for time derivation, the dynamics of x,y, z follows : 

z = -fiiy , Tifii = -2^ PQ E rand (A3a) 
x = -£l y , hfl = e + 2^p a E rand (A3b) ■ 
y = Q x + Sliz (A3c) 

To characterize the random fluctuations in time of E ranc i we model its autocorrelation function by < 
E ra nd(t)E ran d(t + t') > t = -%—W(t' + t c ) — 6(t' — t c )1 where 6(t) stands for the Heaviside step function, r c is the 

characteristic time scale of the fluctuations and u/y/f^ the typical scale of the fluctuating part of the Hamiltonian 
H . This means that E ran< i{t) is drawn at random once every r c and can be considered constant over time intervals 
[nr c , (n + l)r c ], where n is an integer. Within each of these intervals, E ran <i{t) takes the constant value E n . This 
allows to solve exactly the equation for y obtained from Eqs. (A3) : y + (fij^ + ^i,n)V = 0- This yields : 



y(riT c + t) = y(riT c ) cos fl n t H sin Sl n t, (A4) 

where fi„ = ^/fi 2 ^ + fi 2 n with fi ,n and fii ; „ defined as in Eqs. (A3) by setting E ranc [(nT c + 1) = E n . Inserting 
Eq. (A4) into Eq. (A3a) and Eq. (A3b), with the notation X n = X(tit c ) for any quantity X, we get : 

x n+l — X n tt^~ Sn (1 c n) (^-5) 

z n+l — Z n S n 1 ^ V " (1 ~ c ri) (-^6) ) 

ijn+1 = fio^ri+l^ri+l + fil,n+l^n+l (-^7) 

where s n = sinfi„r c , c„ = cosfi„r c . The four equations (A4)-(A7) allow to deduce x, y, z at step (n + 1) provided 
the corresponding quantities are known at step n. Choosing the initial conditions X\,yi, Z\, yields y\ — fi a;i + fii^i 
which allows to initiate the recurrence. Finally, let us note that choosing the initial quantum state as |$i >= 
a\\+ > +yl — |»i| 2 exp(i(pi)\— >, where |+ >, |— > are the eigen states of the TLS, amounts to setting : x\ — 
\a\\yf\ - |ai| 2 cos 991, 2/1 = \a\\*J\ - |ai | 2 sin^i, z\ = \a x f - 1/2. 

Figure 4 shows the dynamics of a TLS defined by A = 1 K, Ao = 0.01 K evolving from the initial state a\ = 
1/2; if i = 7r/2, i.e., from x\ = 0; t/i = \f3j^z\ = —1/4. The random field characteristics were set to u/^/Tc = 0.1 
K and r c = /i/(4e), i.e., r c was chosen four times lower than the Bohr period. Without 'noise', y(t) exhibits the 
well-known regular Bohr oscillations (short-dashed line on Fig. 4). The effect of 'noise' is to deform these ocillations 
(continuous line on Fig. 4) by an amount increasing with time : as a result the periodicity of y(t) gradually disappears. 
This is illustrated in the inset of Fig. 4 showing the exponential decrease in time of the absolute value \C y \ of the 
autocorrelation of y, defined by C y (t) =< Sy(t')Sy(t' +t) > t > /A 2 with Sy(t) = y{t)— <y > and A 2 =< (5y) 2 >. 

Since the value y n depends on the set of values E n drawn for the considered TLS from n = 1, ensemble averaging 
(over many TLS's with the same A, Ao) will lead to a cancelation of y due to the absence of correlations between the 
noise series seen by different TLS's. This cancelation happens on a time scale Ti which should be of the order of the 
one of C y shown in the inset of Fig. 4. This cancelation of y after ensemble averaging amounts to a supplementary 
relaxation term S v /t2 in the Bloch equation describing S y dynamics. 

The dynamics of x(t) (unreported on Fig. 4) is similar to the one of y, yielding a corresponding relaxation term 
S x / T 2- This contrasts totally with the dynamics of z(t), depicted on Fig. 4: provided the amount of noise 8H(t) is 
much smaller than the gap e, z(t) stands very close to its initial value z\, even at large times. In fact small fluctuations 
exist, with an autocorrelation decrease similar to the one of C y , but the key point is that \z(t)/ < z > — 1| <C 1. Hence 
the 'noise' does not yield any supplementary relaxation term in the Bloch equation governing the population dynamics 
S z . 
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FIG. 4: Dynamics of a TLS (A = 1 K, A = 0.01 K) submitted to a random electric field (uj^frl — 0.1 K, r c is the quarter 
of the Bohr period h/ \J A 2 + A 2 ,), z, the quantum mean value of S z , is basically constant (solid line with square symbols), 
i.e., mostly unchanged by the random electric field. On the contrary, y, the mean quantum value of S y , is strongly affected 
by random electric field : the periodic Bohr oscillations (short dashed line) seen in the absence of random electric field, are 
progressively distorded when random electric field is present. Inset: As a result, C y , the normalized autocorrelation function 
of y(t), decreases exponentially with time. 



2. Interaction effects with a finite measuring field E. 



When the measuring E is no longer zero, the whole dynamics should be recalculated, with the supplementary 
dipolar Hamiltonian corresponding to E. However, the fact that the measuring frequency w is much lower than 1/t2 
greatly simplifies the problem. Indeed, if u> were zero, taking into account of E would strictly amount to replace A 
by A + po-E : with this new definition of A, all the previous calculations apply, yielding the same relaxation terms 
in the Bloch equations. We will assume that this holds true for finite u>, due to the fact that for the kHz frequencies 
considered here, the experimental values of T2 ensure u>T2 <C 1, even at the lowest T studied in the body of the paper. 



B. Validity of Bloch equations. 

The three Bloch equations Eqs. (3a)-(3c) are valid in the quasilinear response (34). When the electric field becomes 
strong enough, the relaxation terms form a nondiagonal matrix, e.g. a S z /t x ^ z term might come into play in the 
first Bloch equation, and the corresponding Bloch equations are usually named in the litterature Generalised Bloch 
Equations (G.B.E.). However, up to our knowledge, these generalized relaxation terms have been calculated only in 
the case of transverse fields in the rotating wave approximation |36| . This is at odds with our physical situation : i) 
The transverse field case amounts to A = 0, which, by far, is not the case considered here (remember that, due to the 
1/Ao density of states, for most TLS's one has A > Ao) ; ii) The measuring field E ~ cosuit is an oscillating one, not 
a rotating one ~ expiut and the rotating wave approximation would be valid only close to the resonance ui ~ e/h, a 
condition totally irrealistic here due to the extreme smallness of hoj = 2 x 10~ 7 K. 

However, even if they do not apply in our case, one can use the G.B.E. derived in the rotating wave approximation 
for transverse fields to guess qualitatively what could be the influence of the off-diagonal relaxation terms. Two points 
are worth mentioning : 

i) One can easily check that the GBE still yield qualitatively the quantum saturation phenomenon, even if the 
off-diagonal relaxation terms are responsible for quantitative modifications. In particular, it was shown, in the limit 
of infinite E, that the GBE reduce to the standard Bloch equations with = 2ti and that one gets a vanishing 
susceptibility. 

ii) In the GBE, the off-diagonal relaxation times become infinite (i.e. negligible) when r c — > 0, where t c is the 
correlation time of the random field created, on a given TLS, by its neighbors. In the same spirit p7| , in the GBE, T2 
is affected by a multiplicative factor (l + S 2 r,?) where S = |po-E| /% is the Rabi frequency. The order of magnitude 
of t c in glasses was measured only once by Devaud and Prieur |}8j who found r c ~ 10~ 8 s at T = 70 mK with an 
expected r c ~ 1 /T temperature dependence. The E dependence in the relaxation times can be neglected if Sr c < 1 . 
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Aware of these limits, we guess the standard Bloch equations can give a fair approximation as long as |po-E| does 
not exceed 0.1 — 1 mK at low temperature. As an additional remark, the validity domain of our calculations extends 
as t c decreases. 

To summarize, the GBE do not suppress the quantum saturation phenomenon, on the contrary, they are intended 
to quantitatively account for the various measurable quantities in the saturation regime (such as linewidths, etc.). 
The problem of the strong depression of x' when poE is increased from extremely small values up to 10 -4 — 10~ 3 K 
is thus unavoidable and is at odds with Rogge et aVs experiments 0] which were carried out on various glasses and 
showed absolutely no sign of field induced depression of x'(T < T rev ), despite the fact that p^E was varied from 0.05 
mK to 50 mK : the fact that the domain pqE < 1 mK was experimentally investigated is of special importance since, 
as stated above, in this domain, at least, the Bloch equations used here should be valid. 
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